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Abstract. For the Sun, a tight correlation between various activity measures and oscillation fre- 
quencies is well documented. For other stars, we have abundant data on magnetic activity and its 
changes but not yet on its seismic signature. A prediction of the activity induced frequency changes 
in stars based on scaling the solar relations is presented. This seismic signature of the activity should 
be measurable in the data expected within few years. 
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1. INTRODUCTION 

Only a few years after the discovery of low-degree modes of the solar five minute 
oscillations, Woodard and Noyes (1985) showed using data from the SMM mission 
taken during the declining phase of cycle 21 that the oscillation frequencies decrease 
with decreasing Sun's activity. The correlation between various activity measures and 
the frequencies was definitely confirmed with the data taken during the next activity 
cycle (Libbrecht & Woodard 1990; Bachmann & Brown 1993; Regulo et al. 1994; 
Elsworth et al. 1994). High precision frequency measurements are now available for 
a number of solar-like stars (e.g. Bedding & Kjeldsen 2007). Although, the expected 
seismic signature of the magnetic activity is very weak indeed, its discovery seems to be 
only a matter of time. 

The relative frequency changes in the eleven-year cycle are at most of the order of 
10^^, which is one order of magnitude less than the changes in the irradiance and three 
orders less than in the chromospheric emission measures. We should expect a similar 
situation in distant stars. For many solar-like objects, cyclic variations variations of 
the chromospheric emission are well established (e.g. Baliunas 1998). Nonetheless, an 
effort to measure the much weaker seismic signature is of interest because it reflects the 
changes taking place in deeper layers. 

In expectation of data from measurements, Metcalf et al. (2007, hereinafter M07) 
developed a method for predicting frequency shifts based on scaling the measured p- 
mode frequency variations and changes of the chromospheric activity indices in the Sun. 
Specific forecast was made for the radial modes in the subgiant /3 Hyi. In this paper, I 
outline the principle of the forecasting method adopted in that work and expand on the 
part concerning nonradial modes. 



2. ACTIVITY RELATED CHANGES IN OSCILLATION 

FREQUENCIES 



The starting point in the forecast presented in MOV was the variational expression for 
frequency shifts 

where j = {n,£,m) identifies the mode, is the kernel, which depends on the un- 
perturbed stellar model parameters and eigenvectors of its free oscillations, ^{x). The 
quantity 

Ij = j jVIIP =R^pI (2) 

denotes the mode inertia. The second equality defines the dimensionless mode inertia, 
/. The source, S^, must include the direct feed-back from the magnetic field perturbed 
by an oscillation mode but also the indirect field effects caused by the perturbation of 
the equilibrium density, temperature, and velocity of convection. All these effects may 
be important for frequency shifts at a level of 0. 1 /iHz. Disentangling these effects in the 
solar data is impossible and calculation of their relative contributions to the source is still 
very uncertain. The general form of the kernels for various effects (see e.g. Goldreich 
et al. 1991, Dziembowski & Goode 2004) is complicated but for all of them the leading 
terms are proportional to |div^ 

In their semi-empirical determination of S^{x) for the Sun, MOV adopted the common 
kernel 

j^j{x)^\dHj\^=qjmYn^ (3) 

where D denotes the depth below the photosphere, absorbing the model dependent 
parameters such as gas density or adiabatic exponent in the unknown source rY. Only 
the part of y that is symmetric about the equator affects mode frequencies. Therefore, 
without loss of generality, they could use the following expression 

^W = £^,(D)/'2,(Ai), (4) 

where jl = cos0. It follows from equations ([1]), (|3]), and ([4]) that only the terms up to 
k — £ contribute to the shifts of modes with degree £. 



3. FREQUENCY CHANGES IN THE SUN 

Most of what we know about the changes in the solar oscillation frequencies comes from 
disc-resolved measurements. The data from such measurements are now commonly 
represented in the form 

v„te = £a,-„^.^f(m), (5) 

i=0 

where l3^f{m) are orthogonal polynomials of order i (see e.g. Schou et al. 1994). 
The mean frequencies, uq, yield information on the radial structure of the Sun. The 



higher even-order ayt-coefficients yield separately information on the Sun's asphericity 
described by the individual Legendre polynomial Pikiji) ■ Thus, from measured shifts, 
Aa2i {i = 0, 1, ...), as functions of the mean frequency, we get our constraints on yni^D) 
in the Sun. 

Using equations ^ and (HI) in ([T]), we get 



AV„/m = — — (6) 



where. 



and 



n,lm = j d^d^\Yp\'P2k{ll) = ^^ik{m)Z^ 
f = (-l) 



^(2^-1)!! (2£+l)!! (£-1)! 



This equation with derived from data extending from £ = to over £ = 100, allows 
us to determine the shifts of the low-^ modes much more accurately than if they may 
have been measured directly. 

Data on the frequency shifts show that the source is localized near the Sun's photo- 
sphere. Thus, the following simple form of the depth dependence of the coefficients in 
dH) was assumed in MOV: 

= 1.5 X 10-"A,5(D-D,,,)AiHz2, (7) 

where and D^^k were determined by fitting equation (7) to the data of the a2k coeffi- 
cient shifts, using the relation 

^2k,ne^^kZk,eQn£{Dc,k)^ (8) 

where we denoted 

e„,= 1.5x10-11^. (9) 

The relation ([8]) follows from equations ([5]), ([7]), and Q. The numerical factor is 
arbitrary. M07 determined only the value of Aq, using data for p-modes with n> \ and 
V > 2.5. They showed that the mean frequency shifts, Aao, for over 1500 p-modes are 
well described with our two-parameter fit and that D^q ~ 0.3Mm. There was a small 
systematic departure at low frequencies, indicating that the source extends beneath 0.3 
Mm, but it was regarded to be unimportant. For other a2k of low k, the fit is also good at 
similar values of Dc,k- Thus, data on shifts for each of the low order even a-coefficients 
may be reduced to the single parameters Ak- How these parameters correlate with the 
Mgll activity index is shown in Figure 1. The plotted Ak values were determined on the 
basis of the a2k coefficients calculated by Jesper Schou for the SOHO MDI datcQ. The 
values of AzMgii shown in the bottom panel are from M07 and they were based on the 
NOAA datcll. 



http://quake.Stanford.EDU/ schou/asyml/ 
^ http://www.sec.noaa.gov/ftpdir/sbuv/NOAMgII.dat 
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FIGURE 1. Low order A<. coefficients describing the source of the oscillation frequency changes and 
the changes in the Mgll activity index during activity cycle 23. The coefficient Aq with a la error bar, 
shown in the top panel, describes the mean-frequencies shifts. The next three coefficients are shown in the 
middle panel. The errors bars are similar to those of Aq. 



By comparing the top and lower panels in Figure 1 (similar plots are shown in Figure 3 
of MOV), we see the excellent correlation of the parameter Aq with the changes in the 
Mgll index. Similar correlations exist for the Call and other activity indexes. More 
complicated variability patterns of the higher order Ai^, shown in the mid panel of 
Figure 1, are closely related to the shift of the activity from intermediate to low latitudes. 
We may see that there is a correlation between — Ai and AzMgii but other coefficients 
reach their extreme values at different phases of the cycle. 

Having determined the first four Ay^-coefficients, we may calculate the frequency shifts 
for modes with degrees i <3. Ignoring small changes in the a,- coefficients, related to 
the torsional oscillations, we get from equations ^ and (|7]) 

R ^ 

^Vnlm = T7 E AkQ^f {D,^k) Kkj^, (10) 
k=0 

where R and M are expressed in solar units. We keep these quantities, because we want 
to apply this equation to other stars. 

With known Ayt(?), equation (flOl) allows us to evaluate changes in the mode frequen- 
cies within individual multiplets of low degrees. The plots in Figure 2 show the calcu- 
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FIGURE 2. The changes in the multiplet structure during the activity cycle 23 for two selected mul- 
tiplets calculated from the oq, a\, and 02 SOHO MDl data. The dashed lines refer to modes, which 
correspond to odd values of £+ \m\ and thus are not detectable in the Sun. There is a drastic departure 
from symmetry about the zonal (m = 0) frequency at high activity. 



lated frequency changes during cycle 23 of all components in the two selected multi- 
plets. The plots, which begin at the solar minimum phase, take into account rotational 
frequency shifts, including the linear shift of m x 0.45juHz and the tiny second-order 
effect arising from the centrifugal distortion. At the minimum phase, the multiplets are 
nearly symmetric about the m — mode frequency. The small asymmetry of the £ — 2 
multiplet comes mainly from the 04 coefficient, which arises from the polar field. Even 
the smaller asymmetry of the £ = I multiplet comes solely from the centrifugal dis- 
tortion. At the solar maximum the asymmetry is very large. However, this spectacular 
signature of the magnetic activity is not easy to detect. In fact the asymmetry of an ^ = 2 
multiplet was directly measured by Chaplin et al. (2003), long after the first seismic 
signature of the solar activity has been detected by Woodard and Noyes (1985). 



4. PREDICTIONS FOR OTHER STARS 

The basic assumption in our prediction presented in MOV of the £ — mode frequency 
changes in /3 Hyi was that the relation 

Ao = 22AjMgn5 

which was found for various phases of the solar cycle is valid for other stars. Instead 

of iugih one may use any other global activity index available for the Sun and for a 
considered star. The richest data on global activity are from the Mount Wilson Observa- 



tory - HK Project, which uses Call. Unfortunately, this is a northern hemisphere project, 
which neither includes /3 Hyi, nor any other well studied solar-like pulsators. To apply 
relation (flOl) . we still had to scale Dc o and for that we assumed that it is proportional to 
the pressure scaleheight at the photosphere, that is 

Dcfio^HpOcL^-^^R^-^M \ (11) 

In order to make prediction for nonradial modes, we have to adopt a scaling of the 
parameters at > in equation ([7]). The simplest choice is 



Ak°<^Aooc A/Mgii and Dcj°^Dc,o. (12) 

The first proportionality amounts to assume the same butterfly diagram pattern as in the 
Sun. The second is another extrapolation of the property found in the solar data. Once 
we make these additional assumptions, we may calculate, with the use of equation (flOl) . 
the expected frequency shift of low degree modes for any star for which the shift in iugih 
or any other global activity measure, is known. 

If we had measurements of the individual mode frequencies within multiplets, we 
could test our prediction based on equation (flOl) with single-season observations because 
the same equation describes the momentary v(m) -dependence induced by the momen- 
tary magnetic activity. However, such measurements may be difficult, as we learnt from 
the case of the low-£ solar oscillation. The problem is that the width of the peaks in the 
oscillation spectrum is larger than the frequency splitting between the modes. We thus 
expect that the mean (averaged over multiplet components) frequency shifts will be the 
first detected seismic signature of the magnetic activity in distant stars. 

In order to calculate the mean frequency shifts, we have to weight the individual 
components which depend on the inclination of the rotation axis to the line of sight, Oq. 
The weights should be proportional to the power, which -on average- is proportional to 
\Yp{6o,0)\^. Hence (see equationfTO]), the mean shifts are given by 



2i+lR 

AVnf = > 



2 ^.=0 



AkQne{D,^k)- (13) 



An unknown inclination introduces additional uncertainty to our prediction for nonradial 
modes. 



5. THE CASE OF jS HYDRI 

The G2IV star /3 Hyi was chosen by MOV for predicting radial mode frequency changes. 
Modes with i <2 were measured in the two well separated seasons: 2000.5 (Bedding 
et al. 2001) and 2005.7 (Bedding et al. 2007). There are eight identified p modes with 
£ <2m both data sets. Bedding et al. (2007) find only an insignificant mean frequency 
shift of 0. 1 ± 0.4/iHz between these dates. 

Archival lUE spectra taken between 1978.5 and 1995.5 are available for this star 
allowing to determine the changes in the Mgll index. The new analysis presented in M07 
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FIGURE 3. Predicted oscillation frequency shifts in j3 Hyi between its activity minimum (2004.5) and 
maximum (2010.5). The dates were inferred in MOV on the basis of the i^gn behaviour between 1978.5 
and 1995.5 



has led to an estimated cycle period of 12 y and to a full amplitude of A?'Mgii = 0.15. Both 
numbers are not very different from their solar counterparts. The last observed activity 
minimum took place around 1992.5. In 2000.5 the star was about four years before the 
next estimated minimum and in 2005.7 about a year after. This means that the interval 
is not well suited for testing our predictions but there is no better data available on 
frequencies and activity for any object. 

In order to apply equation (fT3l) for predicting frequency shifts, we need a stellar 
model. The one used by MOV and here has mass, M = 1.1, radius, R= 1 .82, luminosity, 
L = 3.17 (all in solar units), age=6.7 Gy, and the metallicity parameter Z = 0.0138. It 
reproduces the observed mean parameters of the star and the large frequency separation. 
The values of and Dc were calculated according to relations (fTTI) and (fT2l) . adopting 
AzMgii = O.8AzMgii,0 and Dc,q = 0.3 Mm. The numbers calculated for modes identified 
in the /3 Hyi oscillation spectrum by Bedding et al. (2007) are shown in Figure 3. From 
the numbers in this figure, we may see that within a few years we have a chance to detect 
a marginally significant seismic signature of the high activity state of /3 Hyi. However, a 
major improvement in frequency measurements is needed in order to learn more about 
the star's activity changes by means of asteroseismology. Once it is achieved and the 
differences in frequency shifts between modes of different degrees are measured, we 
may get a useful constraint on the longitudinal dependence of the activity changes. 



6. DISCUSSION 



The presented method of calculating an expected seismic signature of stellar magnetic 
activity consists in extrapolating the solar relation between changes in the p-mode 
frequencies and the chromospheric activity indices. The extrapolation lacks a sound 
physical justification so it should be regarded as a very uncertain forecast. 

Numerical predictions of the frequency shifts were given for /3 Hyi, which, like our 
Sun. is not a very active star. The star is brighter and more evolved than the Sun but 
there are no large differences between the two objects. The predicted shifts were found 
similar to those measured in the Sun. They are small and consistent with the null results 
of existing measurements but they may be measurable within the next few years. 

It could be more interesting to test the prediction for a more active object. Unfor- 
tunately, there are none among those stars for which data on solar-like oscillations are 
available. Magnetic activity inhibits precise frequency measurements. Indeed, the anal- 
ysis of the data collected during the solar cycle 23 shows that the widths of the peaks 
in the oscillation spectra increase by some 15% at the high activity phase (Jimenez- 
Reyes et al. 2003). On the other hand, since more active stars are also faster rotators, the 
multiplet structure should be easier to resolve. Thus, I believe such objects should be 
included in the list of asteroseismic targets. The forecast based on the extrapolation of 
solar data to more active stars is more risky. However, regardless whether the forecast 
is right or wrong, its comparison with reality will teach us something interesting about 
stellar activity. 

ACKNOWLEDGMENTS 

This work was supported by the PoUsh MNil grant No. 1 P03D 021 28. 

REFERENCES 

Bachmann K. T., Brown T. M., 1993, ApJ, 41 1, L45 

Baliunas S. L., Donahue R. A., Soon W., Henry G. W., 1998, in Donahue R. A., Bookbinder J. A, eds, 
ASP Conf. Ser. Vol. 85, Cool Stars, Stellar Systems, and the Sun. Astron. Soc. Pac, San Francisco, 
p. 153 

Bedding T. R. et al., 2001, ApJ, 549, LI 05 

Bedding T. R. et al., 2007, ApJ, 663, 1315 

Bedding T. R., Kjeldsen H., 2007, Comm. Asteroseis., 150, 106 

Chaplin W. J. et al., 2003, MNRAS, 343, 343 

Dziembowski W. A., Goode P. R., 2004, ApJ, 600, 464 

Elsworth Y. et al., 1994, ApJ, 434, 801 

Goldreich R, Murray N., WiUette G., Kumar P, 1991, ApJ, 370, 752 
Jimenez-Reyes S. J., Garcia R. A., Jimenez A., Chaplin W. J., 2003, ApJ, 595, 446 
Libbrecht K. G., Woodard M. E, 1990, Nat., 345, 779 

Metcalfe T. S., Dziembowski W. A., Judge P G., Snow M., 2007, MNRAS, 379, L16 (M07) 
Regulo C, Jimenez A., Palle P L., Perez Hernandez R, Roca Cortes T., 1994, ApJ, 434, 384 
Schou J., Christensen-Dalsgaard, J., Thompson M. J., 1994, ApJ, 433, 389 
Woodard M., Noyes R., 1985, Nat., 318, 449 



